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Gauge-spin-space rotation invariant vortices 
in spin-orbit coupled Bose-Einstein condensates 
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We revisit ground states of spinor Bose-Einstein condensates with a Rashba spin-orbit coupling, 
and find that votices show up as a direct consequence of spontaneous symmetry breaking into a 
combined gauge, spin, and space rotation symmetry, which determines the vortex-core spin state 
at the rotating center. For the continuous combined symmetry, the total spin rotation about the 
rotating axis is restricted to 2tt, whereas for the discrete combined symmetry, we further need 2F 
quantum numbers to characterize the total spin rotation for the spin-_F system. For lattice phases 
we find that in the ground state the topological charge for each unit cell vanishes. However, we 
find two types of highly symmetric lattices with a nontrivial topological charge in the spin-^ system 
based on the symmetry classification, and show that they are skyrmion crystals. 

PACS numbers: 67.85.Fg, 03.75.Mn, 05.30.Jp, 67.85. Jk 



C/2 



o 



> 

o 

en 

o 
en 



X 



I. INTRODUCTION 



Vortices are ubiquitous in a rich variety of systems 
from Bose-Einstein condensates (BECs) [ll, 0, superfluid 
helium Q to superconductors [J]. In scalar BECs and 
superfluid helium-4, vortices are quantized, in units of 
a topological invariant called mass circulation ^p v™**^^ • 
dl = 2Trhny/M, where F is a close loop enclosing the vor- 
tex, v™''^'' is the superfluid velocity, and n„ is an integer. 
When internal degrees of freedom are involved, topologi- 
cally stable fractional vortices emerge in the polar phase 
and the cyclic phase of spinor BECs due to the underly- 
ing J.2 [il, [3 and tetrahedral symmetries @ , respectively. 
To classify or predict these and other topological excita- 
tions in spinor BECs such as monopoles, skyrmions and 
knots, we can invoke the homotopy theory Q, because 
the spin and orbital degrees of freedom are decoupled. 

To define topological invariants to characterize quan- 
tized vortices and the quantization of vortices, we usually 
require that the order parameter at the boundary belongs 
to a prescribed order parameter manifold of a certain 
topological nature. However, knowing the boundary of 
the vortex is, in general, insufficient to predict the or- 
der parameter close to the vortex core, because vortices 
havingthe same topological number can have different 
cores [a, Q ■ Therefore, to distinguish vortices with dif- 
ferent types of cores, we should resort to other methods. 
One such method is presented by Salomaa and Volovik 
in Refs. |10l . Ill| . where they applied a symmetry clas- 
sification scheme to distinguish different vortices within 
the same topological class of the superfiuid '^He B-phase 
under an axisymmetric boundary condition, according to 
the axisymmetry and three discrete internal symmetries. 

Here, we show that the symmetry classification scheme 
can also be applied to distinguish vortices in BECs in 
which the spin and orbital degrees of freedom are cou- 
pled. For spin-orbit (SO)coupled spinor BECs [H, [H 
and spinor dipolar BECs |14| - [l6j , the spin and orbital de- 
grees of freedom are coupled due to the single-particle SO 



coupling or the long-range dipolar interaction, and there- 
fore vortices can spontaneously emerge in the ground 
states. Examples are two different vortices found in 
ground states of a spinor dipolar BEC [13, ll3 : with their 
cores filled by the polar phase or the ferromagnetic phase. 
For a SO-coupled spinor BEC, ground states with a single 
vortex or vortex lattices have been found |19l - {25| . 

In this article, we systematically investigate vortices in 
the ground states of a SO-coupled spinor BEC from the 
point of view of spontaneous symmetry breaking (SSB) 
and find that they are direct consquences of SSB into a 
combined gauge, spin, and space rotation symmetry, as 
the system is invariant under the U(l) gauge transforma- 
tion, the combined spin and space rotation and the time 
reversal. The vortcx-core spin state is then determined 
by the spin-gauge symmetry. As long as the spin and or- 
bital degrees of freedom are coupled, possible combined 
symmetries are reduced, because the spin and space must 
be rotated simultaneously to make the Hamiltonian in- 
variant. For lattice phases found in the ground state, 
vortices with different combined symmetries form lat- 
tice structures, and the topological charge for each unit 
cell vanishes. However, in a SO-coupld spin-i BEC, we 
find highly-symmetric skyrmion crystals, which are con- 
structed from superposition of several plane waves. 

Our symmetry classification method differs from that 
in Refs. [ifl[ll| in the following respects. In Refs. |10l . [l]| . 
the authors consider vortices in the superfluid "^He-B 
within the same topological class by assuming an ax- 
isymmetric boundary condition. For SO-coupled BECs, 
we do not take such a boundary constraint. If a vortex 
preserves the continuous combined symmetry, the total 
spin rotation about the rotating axis is 27r. However, if 
the vortex is invariant under the discrete combined sym- 
metry, the symmetry alone is not enough to determine 
the total gauge and spin-state rotations about the rotat- 
ing axis, and we need one additional quantum number to 
count the total gauge rotation and 2F quantum numbers 
to characterize the total spin rotation for the spin-i^ sys- 
tem. Therefore, vortices preserving the continuous and 



discrete combined symmetries are not always fall into the 
same topological class. From numerical results shown 
below, we find vortices preserving the discrete combined 
symmetry and cannot be adiabatically deformed into vor- 
tices which are invariant under the continuous combined 
gauge, spin, and space rotation. 

This paper is organized as follows. Section II describes 
the model Hamiltonian for SO-coupled spinor BECs. 
Section III classifies vortices from their combined symme- 
tries and determines the corresponding vortex-core spin 
state at the rotating center. Section IV discusses different 
types of vortices which are invariant under the continous 
and discrete combined symmetries found in SO-coupled 
spinor BECs. Section V extends the symmetry classifica- 
tion of vortices into spinor BECs. Section VI summarizes 
the main results. 



II. MODEL HAMILTONIAN 

We consider a quasi-two-dimensional spin-i^ BEC with 
N atoms and a Rashba SO coupling. The effective Hamil- 
tonian is "H = Ho + Hint- Here, the single-particle Hamil- 
tonian is given by 



-Ho^ Jdpi^^ 



2A/ 



Vo + jiPxF^ + PyFy) 



V', (1) 



where V' = (V'f, • ■ • , V'-f)^, P = {x,y), and M, Vo = 
Moj^p^ /2, and Fx,y are the atomic mass, the trapping 
potential, and the spin-i^ matrices, respectively. The in- 
teraction part of the Hamiltonian depends on the atomic 
spin and it is assumed to take the same form as that 
described in Ref. Um , where 
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Because the singlet-pairing interaction arises only for the 
spin-2 case, we choose 7 = for the spin-2 and spin-1 
cases. 

The above Hamiltonian possesses the full symmetry 
G = U(l) X S0(2) X r, which describes the U(l) gauge 
transformation, the S0(2) simultaneous spin and space 
rotation about the z-axis, and the time reversal, respec- 
tiveh^ According to the symmetry classification scheme 
[26l429j , we classify ground states by their remaining sym- 
metry described by the group isotropy H [2J|. When 
the state preserves the continuous or discrete combined 
gauge, spin, and space rotations, usually a vortex appears 
at the rotating center, with the vortex core completely 
or partially determined by the corresponding spin-gauge 
symmetry. This is similar to but different from the sit- 
uation discussed in Ref. [26|, where the combined gauge 



and space rotation symmetry plays an important role in 
determining zeros of the superconducting order parame- 
ter of heavy-fermion systems; in the present case, there 
are spin degrees of freedom and vortices can be nonsin- 
gular for SO-coupled BECs. 



III. SYMMETRY CLASSIFICATION OF 
VORTICES 

A. Salomaa and Volovik's method 

In this subsection, we briefly review the symmetry 
classification scheme on vortices in the superfluid ■^Hc- 
B phase by Salomaa and Volovik [10, [lH . The procedure 
of this scheme is as follows: 

(1) Firstly, an axisymmetric boundary condition is as- 
sumed. Far away from the vortex core, the order pa- 
rameter is assumed to be Aaj{p = oo, ip, z) = CRajC^"^'^, 
where p, ip, and z denote cylindrical coordinates and Raj 
describes the relative rotation of the spin (a) and orbit 
(j). By approximately choosing the spin coordinates, 
we can set Raj ~ 6a j- Therefore, the operator describ- 
ing the axisymmetry is given hy Q = L^ + Sz — rnl, 
where L^ = —id /dip -\- Fl,"' is the total operator of the 
^-component orbital angular momentum, including the 
external and internal (Fi,"') orbital rotations, Sz is the 
^-component of the spin operator, and / is the identity 
operator. 

(2) Secondly, if vortices preserve the continuous com- 
bined gauge, spin, and space rotation symmetry, the or- 
der parameter should satisfy QAaj ~ 0. This equation is 
met if Aaj takes the following form: 



Ao., = E ^A.-(/'' ^' ^)A^A,^e^("-^-'')^ 



flL/ 



" int 



(5) 



and 



where AJ and A^ are eigenfuctions of operators F" 

Sz with the correponding eigenvalues ly and fi, respec- 
tively. From Eq. ([5]), we can infer that ii m — fi — v ^ 0, 
C^u{p = 0) = 0. Therefore, the continuous combined 
symmetry reduces the allowed vortex-core state. How- 
ever, when different values of ^ and v all satisfy m — ^—v, 
the vortex core cannot be uniquely determined. In 
Ref. [lOl, Salomaa and Volovik consider three discrete 
internal symmetries Vi , V2 and V^ to further reduce the 
ambiguity of the vortex core, with 

Pi^T'e""", V2^TLy^^Sy,^, ,V3=Vir2, (6) 

where V is the parity transformation p — )■ — /?, Vi denotes 
the combined gauge and the n space rotation about the 
z-axis, and V2 denotes the combined time reversal, the 
TT space rotation {Ly^T^) and tt spin rotation (Sy^T^) about 
the y-axis. 

(3) Thirdly, Salomaa and Volovik also consider vor- 
tices break the axisymmetry under the same axisymmet- 
ric boundary condition as above. To make a vortex break 



axisymmetry, they superpose eigenfunctions of the op- 
erator Q having different eigenvalues. For example, if 
the order parameter takes the form A^j = Aa,j{Q ~ 
0) + A„j(g = 2), where 

and C^i/(Q 7^ 0,p = 00) = 0, the corresponding vor- 
tex is invariant under e''^'^. However, as the boundary 
condition is fixed, the phase winding number for each 
component is fixed to be m — ^ — :/ and this vortex is 
topologically equivalent to the one preserving the contin- 
uous combined symmetry. 




FIG. 1: (color online). Schematic spin textures for two differ- 
ent vortices in a spin-i system preserving the same discrete 
combined symmetry described by e"'^'^"Cz(27r/3), where the 
total spin rotation about the rotating axis at O is 2-k for (a) 
and — 47r for (b). 



B. Our method from a SSB point of view 

In this subsection, we systematically classify vortices 
in SO-coupled spinor BECs from a SSB point of view. 
When the space rotation is involved in the SSB, a vortex 
might emerge. In general, to classify the symmetry of 
vortices, we should consider different subgroups of the 
full symmetry group G of the Hamiltonian. Due to the 
SO coupling, we need to simultaneously rotate the same 
angle of the spin and the space, which therefore simplifies 
the procedure to find subgroups of G. 

Firstly, we neglect the time-reversal symmetry. Then, 
there are two different remining symmetries, that is (I) 
continuous symmetry and (2) discrete symmetry. The 
corresponding generator is given by e^'^Cz{0)^ where Cz{0) 
refers to the combined 9 spin-space rotation about the z- 
axis. For the states with continuous symmetry we have 
e e [0,27r) and = k6' with k = 0,±1,±2,... (k = 
±1/2, ±3/2, . . .) for the integer (half-integer) spin case, 
whereas for the states with discrete symmetry we have 
9 = 2-K /n (n e Z) and (f> = k9 with k = 0, . . . , n — 1 
(k = 1/2, . . . ,n — 1/2) for the integer (half-integer) spin 
case. 

To preserve the combined gauge, spin, and space rota- 
tion symmetries, the vortex-core spin state \tJj) located at 
the rotating center must satisfy the following spin-gauge 
symmetry: 



iF,d 



^nv') = w- 



(8) 



We first focus on the continuous symmetry. In con- 
strast to the superfluid hclium-3, where the continuous 
combined symmetry is not usually sufficient to uniquely 
determine the vortex core, the order parameter of a SO- 
coupled BEG is described by 2F+1 complex numbers and 
we can uniquely determine the state of the vortex core 
as (1) nonsingular with \tp) ~ \1\If = i^) ii —F < k < F 
and (2) singular with the vacuum core if \k\ > F. For 
instance, in the F — 1 case, two different nonsingular 
vortices correspond to the vortex core filled by the polar 
and the ferromagnetic phases. 



Next, we turn to the discrete symmetry case. Differ- 
ent from the continuous symmetry case, whether we can 
uniquely determine the core state depends on the values 
of K, n and F. For instance, when k = 0, there is only 
one type of nonsingular vortex core with \ip) = \Mp = 0) 
for the integer spin case, whereas when k = 1, n = 3, 
and i^ = 2, we can find that an arbitrary superposition 
of \Mf = 1) and \Mf = -2) satisfies Eq. dH), which 
means that we cannot uniquely determine the nonsingu- 
lar vortex-core state. Table U summarizes the relation 
between the combined symmetry and the vorte-core spin 
state. 

No matter whether we can uniquely determine the 
vortex-core spin state for the discrete combined symme- 
try, there are still ambiguities for the total gauge and 
spin rotations around the vortex core, which is different 
from the case of continuous combined symmetry, where 
the total spin rotation about the rotating axis is 27r and 
the total gauge rotation is 2k7t. Figure [T] schematically 
illustrates two different vortices in a spin-i system which 
preserve the same discrete combined symmetry described 
by e^'^'^^/'^Cz{2'K/3). From this form, we can infer that 
the allowed total spin rotation along the loop PqPi is 
27r/3 + 27rNs (Ns € Z) with (a) N^ = and (b) N^ = -1. 
A similar ambiguity exists for the gauge rotation. There- 
fore, to remove the ambiguity, we can define quantum 
numbers to characterize vortices preserving the same dis- 
crete combined symmetry. 

For the spin-i case, two quantum numbers Ng and 
Ns are needed to remove the ambiguity in the gauge and 
spin rotations, respectively. If the vortice is invaraint un- 
der e*''^'^/"Cz(27r/n), the total gauge rotation and spin 
rotation about the rotating axis are 2kt: + 27rn[l\lc, and 
2tt + 27rnNs, respectively, where Ng — mod(Ns, 2)/2 G Z, 
because if we rotate the half-integer spin for 27r, we need 
to perform a n gauge rotation to make the state invari- 
ant. To determine two quantum numbers from numerical 
results, we can focus on the phase winding number for 



TABLE I: Combined symmetry, vortex generator, and vortex-core state. For the discrete combined symmetry, we consider 
n > 3 for the spin-i and spin-1 cases, and n = 3, 4 for the spin-2 case. For the case of ti = 2, we cannot predict whether there 
is a vortex. When no spin state at rotating center can satisfy the corresponding spin-gauge symmetry of Eq. ((Sjl, singular 
vortices should appear with the vortex core filled by the vacuum/normal state. 



combined symmetry 


vortex generator 


vortex-core state 


continous 


€"""^{6) withSG [0,27r) 


\Mf — n) if h = —F, . . . , F 
vacuum/normal state if |k > F 


discrete (-F = 5) 


e'"''"C.(27r/n) 
e^'^'2^/"C,{2TT/n) with k = 3/2, . . . , n - 3/2 


spin-up 

spin-down 

vacuum/normal state 


discrete (F — 1) 


e'2-/"c,(27r/n) 

C.(27r/n) 

e'("-i)2-/"C,(27r/n) 

gi«2^/n(j^(27r/n) with k = 2, . . . , n - 2 


\AIf — 1) (ferromagnetic) 

\Mf = 0) (polar) 

\AIf = — 1) (ferromagnetic) 

vacuum/normal state 


discrete {F = 2) 


C,(27r/n) 
e'2-/3c,(27r/3) 
e'*"/='C,(27r/3) 

e'3"/2^.(7r/2) 
e"C,(7r/2) and T 


\Mf = 0) (uniaxial nematic) 

arbitrary superposition of \Mf = 1) and \Mf = —2) 

arbitrary superposition of \Mf = 2) and A/f = —1) 

\Mf — 1) (ferromagnetic) 

\Mf = — 1) (ferromagnetic) 

(-g-i,^/2^Q^Q^Q^gi0/2-)T ^j^j^ ^ arbitrary (biaxial nematic) 



the spin component \Mf) defined as 

1 f Vm^^^Mf - {^i'*M^)i^M„ 



WMf 



27r 



2|VAf^Pi 



dl (9) 



where F is a loop around the vortex core, as we have 



wmf = k + 7ilMg — (1 + nMs)Mf 



(10) 



If his 7^ 0, we find a topologically different vortex in com- 
parison with the vortex which preserves the continuous 
combined rotation symmetry. 

Different from the spin-^ case, for the integer spin-i^ 
case, we need 2F + 1 quantum numbers: one quantum 
number Ng for describing the total gauge rotation and 2F 
quantum numbers for denoting the total spin rotation, 
because we need 2F vectors for each spin state in the 
Majorana representation [2^, |3i|. Equivalently, we can 

define another 2F quantum numbers Ns ^ for removing 
the ambiguity of the total spin rotation by evaluating 
phase winding numbers for nonzero spin components. If 
the vortex is invariant under e*"^'^'"C2(27r/n), we have 

wmf =K + n\!^g- (1 + jiNi*^^^") Mp, if Mp + 0, 

it'o = K + nh\g, (11) 



and Wf^"^ X Mf 

|(Mf) 



e Z. If at least one 



where N^ G 

quantum number Ns'"^'' is nonzero, the phase winding 
numbers for all spin components cannot be written in a 
compact form as wmf = —Mp -I- Nc, with Nc S Z, which 
is the salient difference compared with vortices preserving 
the continuous combined rotation symmetry. 

Taking into account the time-reversal symmetry, sim- 
ilar arguments can be applied. In special cases, the al- 
lowed vortex-core spin state is further restricted due to 



the time-reversal symmetry. For instance, for the spin-2 
case, if a vortex preserves the discrete combined symme- 
try described by e^'^Cz{T^/2), the vortex-core spin state 
can takes an arbitrary superposition of \Mp = 2) and 
\A'Ip = —2). If we require the vortex be invariant under 
the time reversal, the allowed vortex-core spin states are 
fmther reduced, as listed in Table HI 



IV. VORTICES IN SPIN-ORBIT-COUPLED 
BECS 

In this section, we show how to understand different 
vortices found in the ground states of SO-coupled BECs 
from the standpoint of combined symmetries. 



A. Continuous combined symmetry 

We first focus on the vortices invariant under the con- 
tinuous combined gauge, spin, and space rotation symme- 
tries. In Refs. |2J], approximate low-lying single-particle 
eigenstates of the Hamiltonian "Ho with strong SO cou- 
pling are presented in momentum space as 

V'f..™(k) oc k'-^'^e-^'^'-^'^'/^Hr^ik' -v') 

xe™^''C-i=^(<^fc), (12) 

where k is the wave vector, Lp^ = arg(fc2: + iky)^ a±_ = 
y^h/Muj± , v' = v/uj±a±, k' ~ ka±, Hn are Hermite 
polynomials and 

C_i(^fe)-(l,-e^^'^-)^/A 

C-i(¥>fc) - (e-^'^^-V2,e*^'=)^/2, 
C-2(v5fe) = (e-2''^^-2e^*^^ V6,-2e*^^e2*^'=)^/4.(13) 
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FIG. 2: (color online). Ground-state spin textures for three different lattice phases showing triangular (a) and kagaome (e, g) 
structures for a spin-orbit-coupled pseudo-spin-i BEC, with v' = 15 and (a) a = 0.5hiu±/N and /3 = —0.1a, (e) a = 0.3hLu±/N 
and /3 = 0.1a, and (g) a = OAhuj±/N and /3 — 0.1a. The corresponding spin textures for the variational order parameters 
constructed by a superposition of several plane waves given in Eq. (|15|l are shown in (b), (f), and (h), respectively. Here 
the spin polarizations along the x-y plane are denoted by arrows and its z-component (a.) is shown according to the color 
gauge on the right bottom. The phase distributions of (c) spin-up and (d) spin-down components for the ground state of 
the triangular-lattice phase in (a) are also shown according the color gauge on the left bottom. The size of each figure is 
[—ax/2, ax/2] X [—ax/2, ax/2]. A unit cell of the order parameter is indicated by a solid hexagon. Singular vortices indicated 
by open circles form a triangular lattice for all cases. Two nonsingular vortices are indicated by filled circles with blue (^3) 
and red (^4) colors. 



By the Fourier transformation, we obtain the correspond- 
ing order parameters in real space, e.g. for the spin- 2 
case, it takes the following form: 



Jo 



i'"+V„+i(Me'("+'^^ j "^^'^^^^ 



where (p — aTg(x+iy) and Jm is the Bcssel function of the 
first kind. All states given in Eq. ([T^ possess the conti- 
nous gauge, spin, and space rotation symmetries. Some 
of them can be stabilized with weak interactions. For 
instance, two different vortices are found in the pseudo 
spin- 2 case, and their order parameters take the form 
given in Eq. ((Ti)) with n = and m = —2,-1,0,1. 
They are invariant under e^'^^Cz{9) with k = ±1/2 (see 
Fig. 3(d) of Ref. [H and Fig. l(a)of Ref. [H) and 
K ^ ±3/2 (see Fig. 3(b) of Ref. ^ and Fig. 1(b) 
of Ref. [23|)- From the combined symmetry, we know 
that the vortex-core state at the rotating center should 
be invariant under the spin-gauge symmetry g^K^g-^^^^, 
Therefore, we determine the vortex core as that sum- 
marzied in Table ID 



We would like to point out that such vortices with 
continous combined symmetry also appear in a spinor 
dipolar BEC where the effective SO coupling arises from 
long-range magnetic dipole interaction. For instance, two 
different ground states preserving continous combined 
symmetry are found in the ground state of the spin-1 
^^Rb condensate [13, ll3 ■ One called the chiral spin vor- 
tex phase is invarint under e^^Cz{d), which results in a 
vortex with the vortex core filled by the ferromagnetic 
phase. The other called the polar-core vortex phase pre- 
serves Cz{d), which induces a polar-core vortex because 



-iF^t 



\Mf = 0) = \AIf = 0)- 



B. Discrete combined symmetry 

According to the symmetry classification scheme [241 . 
there are states preserving discrete rotation symme- 
tries. Some of them spontaneously emerge as the ground 
states with individual spin-component density distribu- 
tions showing lattice structures [20|-[24|. In the strong 
SO-coupling limit, the order parameter can be well ap- 
proximated by a superposition of several plane waves 



given in Eq. (6) of Ref. [2J] as 



3=0 



(15) 



where hj = fcg(cos(j27r/[l\l),sin(j27r/[l\l)) and kg = mv/h. 



1. spm-1/2 

Firstly, we focus on the pseudo spin-^ case, where 
three different types of lattice phases are found in the 
ground state. Figure [2] illustrates their spin textures 
which are obtained by numerically solving the coupled 
Gross-Pitaevskii equations with v' = 15 and (a) a = 
0.5hLo±/N and /3 = -0.1a, (e) a = 0.3huj±/N and 
/3 = 0.1a, and (g) a = OAHuj±/N and /3 = 0.1a, 
and compares them with those of the variational order 
parameters described in Eq. (|15|) with parameters (b) 
N = 3 and ((/>o, </>!, 02) = (tt/S, tt, 57r/3), (f) IN = 6 and 
('/'o,'/'i,02,03,04,'/>5) = (7r/6,7r,57r/6,-7r/3,-7r/2,7r/3), 
and (h) N = 6 and (f>j = {j + l)7r/3 with j = 0, . . . , 5. 
The validity of the ansatz in Eq. (|15p is numerically con- 
firmed. In the following, we discuss the symmetry of the 
vortices based on the ansatz. Strictly, the symmetry we 
discuss below can only be approximately satisfied. If we 
increase the strength of SO coupling, the agreement be- 
tween the ansatz and numerically exact result becomes 
better. 

The triangular-lattice phase shown in Figs. [2l^a-d) is 
found to be composed of half and anti-half skyrmions 
(enclosed in the triangle A0A1A2 and the one nearby), 
which appear alternatively and are connected by singular 
vortices (with one locating at Aq). Based on the symme- 
try analysis in Sec. IIIH we find that for the case of n = 3 
there are three different types of symmetries generated 
by e*'^27r/3(-^(27r/3) with k = 1/2,3/2,5/2, respectively. 
As presented in Ref. [2J|, they are connected by shifting 
the lattice as shown in Figs.Hfa-d), where three different 
vortices preserving different combined discrete symme- 
tries are denoted by Aq, A3 and A4, filled by the vacuum, 
spin-up, and spin-down states, respectively. To confirm 
this statement, we focus on the order parameter 



\/3y ^ 



^ = e*'^/3g»fc,x^_^(Q^_^gZ^g^fe,(-f + ^)^^^(27r/3) 



_|_g''57r/3gifcg(- 



•3» 



C-i(M3), 



(16) 



which is shown in Fig. ^h). It is invaraint under 
e"Cz(27r/3). Therefore, vortices at the rotating center 
denoted by Aq is singular. Shifting the lattice center to 
the point ^3 by coordinate transformation {x-|-27r/3fcg — >■ 
X, y — 27r/3\/3fcg — >■ y}, the order parameter is trans- 
formed into 

^ = e*^''/3e*'=''^C-i(0) + e'^''/^e'''s(-f+4^)^_i(27r/3) 



which is invariant under e*'^/'^Cz(27r/3). We thus find 
that the vortex core is filled by the spin-up state. Similar 
arguments can be applied to vortices located at Ai, A2 
and A4. 

In constrast to the triangular-lattice phase, two 
kagome- lattice phases shown in Figs.[2l(e,f) and (g,h) can- 
not be recast into half- and anti-half-skyrmion lattices. 
They possess combined discrete symmetries with symme- 
try groups generated by {e*(2''/3+'^/3)C^(27r/3), rC^(7r)} 
and {e*(''/3+7r/6)^^(^/3)| [3^^ respectively. Therefore, 
we expect that the vortex in the rotating center should 
be filled by the vacuum state, as denoted by Bq and Co, 
respectively. Besides, there are two more different singu- 
lar vortices denoted by Bi{Ci) and i?2(C2), originating 
from different symmetries, which can be determined by 
shifting the lattice center as that applied in discussing 
vortices in the triangular-lattice phase. 

Refocusing on spin textures shown in Fig. [21 we do 
see that for vortices at Aq, Bq, Bi, B2, Ci and C2, the 
total spin rotation around the vortex core is — 47r, dif- 
ferent from 2tt. The corresponding quantized number 
Ng is —1. Such vortices are unique to the discrete com- 
bined rotation symmetry. By examining the phase wind- 
ing number for each spin component, we can determine 
another quantum number N^ by solving Eq. ((TU]). To 
show this explicitly, we take the vortex at Aq preserv- 
ing e"Cz(27r/3) as an example. From Figs. HJc-d), we 
infer that the corresponding phase winding numbers for 
two spin components are wi = 1 and w_i = — 1. Using 
Eq. (Uni) with Mf = 5, we get 



^+3N,-(1-3)Q)=1. 



(18) 



Therefore, we find N„ = — ^ 



2. spin-1 & spin-2 

Next, we focus on the case of integer spin, where we 
find triangular and square-lattice phases. They can all 
be well described by superposition of three or four plane 
waves as described in Ref. [2J| and in Eq. (|15p with INJ = 3 
and 4. Similar to the spin-^ case, the combined sym- 
metry is approximately satisfied for numerically derived 
ground states. In the following, we only use their cor- 
responding variational order parameters in Eq. (|15p to 
discuss the symmetry. Figure 2 shows the corresponding 
spin textures in a unit cell of five lattice phases found in 
SO-coupled spin-1 and spin-2 BECs. For the triangular- 
lattice phase, we can choose arbitrary values of <f>j in 



Eq. (|15p . as they are connected by shifting the lattice 
center. For the square-lattice phase in Figs. [3l^c-e), one 
representative state in Eq. ([15]) can be chosen as (c, d) 
4'i ~ J'''/2 and (e) (f)j = with j = 0, . . . , 3. 

Similar to the pseudo spin-2 case, three different vor- 
tices appear simultaneously in a triangular-lattice phase 
due to three different combined discrete symmetries gen- 
erated by e"*^^/^Cz(27r/3) with k = 0,1,2, respectively. 
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FIG. 3; (color online). Spin textures of the triangular-lattice phase (a, b) and the square-lattice phase (c, d, e) of spin-orbit- 
coupled spin-1 (a, c) and spin-2 (b, d, e) BECs. Each subfigure shows spatial variation of the corresponding order parameter in 
a unit cell visualized by plotting '^^ ^rnYp'{0, ip) with its phase displayed according to the right color gauge, where Yp ^ is the 
spherical harmonics of rank F and ^m is the m-th component of the spin wave function. For the triangular-lattice phase, there 
are vortices with cores at the vertices and the center of the triangles. For the square-lattice phase, vortex cores are located at 
the center of the squares. One representative vortex for each type of symmetry described by e^'^'^'^'^Cz{2TT/3) [e"'''"Cz(7r/2)] is 
enclosed in the circle 'o' and denoted by Pk [Qk]- 



They are connected by shifting the lattice center [2J]. 
By directly solving Eq. (j8|), we uniquely determine the 
vortex-core spin state for the spin-1 case. This breaks 
down for the spin-2 case on two specific symmetries 
of e'2'^/3^^(27r/3) and e'^''/^C,{2n/3), as there are in- 
finite vortex-core spin states satisfying the same com- 
bined symmetry. For instance, the spin state ci\Mf = 
1) -I- C2\Mp = —2) with arbitrary coefficients ci and C2 
is invariant under e'^'^/'^e"'^^^'^/^, and the one which is 
actually chosen is determined by the Hamiltonian. In 
Figs. [Sl^a) and (b), we indicate by circles three different 
types of vortices which preserve different symmetries. 

For the case of n = 4, as discussed in Sec. IIIIBi there 
are four different vortices preserving the combined dis- 
crete symmetries generated by e*"'^/^C^(7r/2) with k = 
0, 1,2,3. They are found in two different square-lattice 
phases with one breaking time-reversal symmetry as illus- 
trated in Figs, ^c) and (d) and the other one preserving 
time-reversal symmetry as illustarated in Fig. ^e) . As 
shown in Ref. [2J] , in a time- reversal symmetry breaking 
square-lattice phase, the state can preserve two different 
combined symmetries of e^'^^'^Czij^ j'i) and e*'^'^/^Cz(7r/2) 
by choosing the proper center of the lattice. This re- 
sults in two different vortices with the corresponding 
vortex-core spin state preserving the spin-gauge symme- 
try e»'^/2e-»^-'^/2 and (^'i^l'^e-^^-'^l'^ . Solving Eq. ([5]), 
we can uniquely determine the vortex-core spin state, as 
demonstrated in Figs. [DJc) and (d). The other type of 
square-lattice phase, which preserves time-reversal sym- 
metry and Cz(7r/2) or e^'^Czij^ jT) relying on the choice of 



the lattice center, has only been found in the spin-2 case. 
These symmetries induce two different nonsingular vor- 
tices with vortex cores determined from Eq. ([8]) and the 
time-reversal symmetry requirement. These two vortices 
are shown in Fig. [3l[e). 

Reconsidering vortices shown in Fig. [3l we find that 
there are new vortices which are unique to the discrete 
combined symmetry with at least one nonzero quantum 

number Wg ^ related to the total spin rotation. To 
demonstrate the searching procedure explicitly, we take 
the triangular-lattice phase of spin-1 BECs as an exam- 
ple. The ground-state phase distributions for three spin 
components are shown in Figs. |4l^a-c). To examine the 
vortex at P\ preserving e*^'^/^Cz(27r/3), we use Eq. (fTTj) 
to derive the following three equations: 



1 + 3N„ 



l + 3Ng- (l + 3Ni^') =0, 



1 + 3N 



1 



-3Ng = 1. 



(19) 



(1) M(-l)^ 



Solving these equations, we find (INg,INs ,N 
(0,0,-1). Therefore, this vortex is topologically differ- 
ent from vortices preserving the continuous combined 
symmetry. Applying the same method to vortices at 
Po and Pi preserving Cz(27r/3) and e^'^''l'^Cz{2'K jZ), re- 
spectively, we find their three quantum numbers as 
(Ng,Ni'\N,^'^) = (0,0,0) and (-1,-1,0). 

As discussed in Sec. IIIIB| the reason why we need 
IF quantized numbers to describe total spin rotation is 
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FIG. 4: (color online). Phase distributions for individual spin 
components {Aip = 1, 0, —1) (a-c) and the corresponding spin 
textures for two unit vectors using the Majorana representa- 
tion (d-e) of the ground state showing triangular structure 
of a spin-orbit-coupled pseudo-spin-1 BEG with v' — 15, 
a — 0.5hLjj±/N and /3 = 0.2a, A unit cell of the order pa- 
rameter is indicated by a solid hexagon. Vortices denoted by 
Po, Pi and P2 are the same as those in Fig. [31 a). 



that to denote each spin-F state, we need 2F unit vec- 
tors, which can be obtained by using the Majorana rep- 
resentation [23,[3l[, where 2F unit vectors (71^,71^,71^) = 
(sin 6' cos (y9, sin sin <p, cos 6*) correspond to the 2F solu- 
tions of the equation 



h^ 



h2Fr = 0, 



(20) 



where ^ 



'tan(6'/2) and h = ip*p_J^/{2F ~ i)W. 



Figures IH^d-e) show the corresponding spin textures of 
two unit vectors for the ground state in Figs. lUJa-c), 
where we classify two unit vectors numerically by com- 
paring values of their z-componcnt n^. We can infer that 
both two vectors rotate around the rotating axis at Pq for 
2tt, whereas for vortices at Pi and P2, the total rotation 
for one vector is —An. This is consistent with the above 
arguments that the vortex at Pq is topologically equiva- 
lent to vortices which preserve the continuous combined 
symmetry, whereas vortices at Pi and P2 arc topologi- 
cally different. 



C. Lattice translation symmetry 

For lattice phases, vortices form lattice structures. Be- 
sides preserving the discrete combined rotation symme- 
try, their order parameters are approximately invariant 
under lattice translation. For such lattice structures we 
can define a topological charge for each unit cell which is 
analogous to the Chern number: 

Ch = // (V X v--^^),da;rfy, (21) 

where the superfluid velocity v"'^'^^ is given by 



2M|V'|^ 



:[V^*VV'-(VV^)*V^]. 



(22) 



For the integer spin cases, the topological charge for dif- 
ferent lattice phases in the ground state vanishes, as 
there are only nonsingular vortices and the order pa- 
rameter changes smoothly over the entire manifold. For 
the pseudo spin-^ case, we have confirmed that all three 
different lattice phases have zero topological charges, al- 
though there are singular vortices. Based on the Mermin- 
Ho relation ^3 



(V X v'"''"") 



AttM 



e^,rS ■ {d,S X drS), (23) 



where S*^ = 



fff^) and {ii,v,t) = x,y,z, we find that the 
total skyrmion number Qskx for each unit cell is zero, 
where 



Qskx = ^ // S • (a,S X dyS)dxdy. (24) 

Therefore, the lattice phases found in the ground state 
of a SO-coupled BEC are quite different from skyrmion 
crystals (SkXs) with a nonzero topolo gica l charge for 
each unit cell found in chiral magnets 1331 . where the 
Dzyaloshinskii-Moriya (DM) interaction [3J] plays a cru- 
cial role. However, we find that for the spin-i case, sev- 
eral lattice phases, which we predict from the symmetry 
classification [2J] but have not been found in the ground 
states, show similar structures as SkXs. Figure [S] shows 
two different SkXs, where the corresponding order pa- 
rameters are given by Eq. (|15p with (a-c) N = 4 and 
(l)j = J7r/2 (j = 0, . . . , 3) and (d-f) N = 6 and 4>j = j27r/3 
(j = 0, . . . , 5). From the spin textures in Fig. [5jc) and 
(f ) , we can infer that both SkXs have the same topologi- 
cal charge Qskx = —2. Here, the unit cell is determined 
from the periodicity of the order parameter. By per- 
forming time reversal, we obtain two other SkXs with 
<9skx = 2. The square shape of the SkX in Fig. [SJa-c) 
has also been found as a stationary state in a system 
without SO coupling |35| . 



V. VORTICES IN SPINOR BEGS 

In this section, we generalize the symmetry classifica- 
tion of vortices to a system in which the spin and or- 
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FIG. 5: (color online). Two different skyrmion crystals predicted from the symmetry classfication for the spin-i case [2J|, where 
their spin-up, spin-down component density distributions, and spin textures are shown from left to right. The corresponding 
order parameters are given by Eq. psp with (a-c) N = 4 and cfij = JTr/2 {j — 0, . . . , 3) and (d-f) N = 6 and 4>j — j2tt/3 
{j — 0, ... ,5). A unit cell is indicated by (c) a solid square and (f) a solid hexagon, both with Q = —2. 



bital degrees of freedom are decoupled. Examples in- 
clude spinor BECs and rotating spinor BECs, where the 
corresponding Hamiltonian is invariant under the U(l) 
gauge transformation, S0(3) spin rotation, and S0(2) 
space rotation about the z-axis. For spinor BECs, we 
can choose their Hamiltonian as a special case in Sec. HIl 
with V — 0. Furthermore, by adding a rotating term 
'Hrot = —ft J dp ijj^Lzipi where fl and L^ refer to the 
rotating frequency and the projected orbital angular mo- 
mentum, we obtain the Hamitlonian to describe a rotat- 
ing spinor BEC. 

As long as the spin and orbital degrees of freedom are 
decoupled, when spontaneous symmetry breaking occurs, 
there are other types of vortices which are invariant under 
the combined symmetry described by e*'^e^*^^-7?.^('i?), 
where 0, 9, and '& refer to the gauge, spin, and space 
rotation angles, respectively, and TZz{'d) describes the 9 
space rotation about the z-axis, because 9 and d can 
take different values. As there are numerous vortices 
from the viewpoint of the spontaneous symmetry break- 
ing, in the following, we focus on a special vortex with 
fixed boundary conditions and predict their possible vor- 
tex cores based on the combined symmetry in a manner 
similar to the classification scheme of vortices in super- 
fluid ^He-B [13, El- 

For instance, we take a fractional 1/2-1/4 vortex in 
the biaxial nematic phase of spin-2 BECs discussed in 
Ref. [31, where 1/2 and 1/4 refer to the tt gauge trans- 
formation and the 7r/2 spin rotation around the vortex, 
respectively. Therefore, away from the vortex core, the 
axisymmetric order parameter can be chosen as ■0(p — >■ 



oo) 

-05 



g«W2e-'¥>F./4^^ = 7pl(l,0,0,0,e*v')^/v^, where 

' \fPs{^,^T^T^A)^ /^/^ is the order parameter for 
the biaxial nematic phase, and ip is the azimuthal angle. 

Under this boundary condition, the continuous com- 
bined symmetry which can be satisfied by the vortex 
is described by operators e^'^^e~^^^'TZz{4:9) with 9 g 
[0,7r/2], which implies that the tt gauge transforma- 
tion and the 7r/2 spin rotation are uniformly distributed 
along the azimuthal direction. When 9 = 7r/2, we get 
gj7rg-.»_F,7r/2-^^(27r)|V') = e^e-*-^-'^/^ | V') = IV'), which 
results in a strong constraint for the order parameter 
of the fractional 1/2-1/4 vortex satisfying the continu- 
ous combined rotation symmetry, because ip should be 
a superposition of \Mp = 2) and \Mp = —2) spin 
states everywhere. Meanwhile, we can determine the 
nonsingular vortex-core spin state at the rotating cen- 
ter from the spin-gauge symmetry, and find that it is a 
ferromagnetic phase with order parameter ip{p = 0) = 
Y^(l, 0, 0,0, 0)"^. This type of fractional vortices filled 
by the ferromagnetic core has already been numerically 
found in Ref. [9[. 

Besides the continuous combined symmetry, we can 
also search for a fractional 1/2-1/4 vortex with a discrete 
combined symmetry, which is described by the opera- 
tor e^'^e~^^^'TZz{2Tr/n). Here, we focus on the discrete 



inip ^ — in9Fz 



t 



2F-I-1- 



symmetries with a constraint (A) e™'^e 
Therefore, the order parameter is no longer fixed into a 
superposition of two spin states. To be consistent with 
the boundary condition, there is another constraint (B) 



gi0g-zeF, 



(1,0,0,0,0)^ 



jW»e-*'^-'^/2"(l,0,0,0,0)^. 



By solving these two equations (A) and (B), we can pre- 
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diet possible types of fractional 1/2-1/4 vortiees distin- 
guished by the discrete combined gauge, spin, and space 
rotation symmetry. One solution is that ~ 47r/3, 
9 = 27r/3 and n = 3. Furthermore, we determine the 
nonsingular vortex-core spin state at the rotating cen- 
ter to be an arbitrary superposition of \Mf = 2) and 
\Mf = —2) spin states. The coefficients of the superpo- 
sition depend on the Hamiltoinan. In Ref. [9| , a specific 
vortex with a cyclic core is discussed. 

We have predicted possible vortices and the corre- 
sponding vortex cores under the axisymmetric boundary 
condition for a fractional 1/2-1/4 vortex. Such symme- 
try classification is also applicable to understand other 
vortices with axisymmetric boundary conditions. 



VI. CONCLUSIONS 

We have systematically studied the structure of ground 
states of SO-coupled spinor BECs, including pseudo spin- 
2, spin-1, and spin-2 systems. In constrast to quantized 
vortices in spinor BECs which are distinguished by top- 
logically invariant quantized numbers, vortices in SO- 
coupled BECs spontaneously emerge in the groud states, 
and are classified by the combined gauge, spin, and space 
rotation symmetry and even time reversal. In addition, 
the symmetry classification not only classifies vortices 
but also determines the vortex-core spin state. We have 
studied vortices with both continuous and discrete sym- 
metries. For vortices preserving the continuous combined 
symmetry, the total spin rotation around the vortex core 
is fixed to be 27r, whereas for vortices invariant under the 



discrete combined symmetry, we need a further quantum 
number Ns for the spin-i case and 2F quantum num- 
bers hli "' with Mp — ±1, . . . , ±F for the spin-F case 
to specify the total spin rotations. From numerical re- 
sults, we do find new vortices with nontrivial quantized 
numbers (N^ ^ 0, N^ '' ^ 0). For lattice phases, where 
the order parameter preserves lattice translation symme- 
try, we define a Chern-number-like topological charge for 
each unit cell. Unfortunately, all lattice phases found in 
the ground state have zero topological charge. However, 
we find two types of skyrmion crystals in the BEC phase 
of the spin-i system. They are highly symmetric states 
predicted from the symmetry classification, and may be 
experimentally realized as stationary or excitated states. 
Finally, we also show that the symmetry classification of 
vortices is also applicable to other systems such as spinor 
dipolar BECs, spinor BECs and rotating spinor BECs. 
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